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ABSTRACT 

The coalescence of massive black hole (MBH) binaries following galaxy mergers is one of the main 
sources of low-frequency gravitational radiation. A higher-order relativistic phenomenon, the recoil as a 
result of the non-zero net linear momentum carried away by gravitational waves, may have interesting 
consequences for the demography of MBHs at the centers of galaxies. We study the dynamics of recoiling 
MBHs and its observational consequences. The "gravitational rocket" may: i) deplete MBHs from late- 
type spirals, dwarf galaxies, and stellar clusters; ii) produce off-nuclear quasars, including unusual radio 
morphologies during the recoil of a radio-loud source; and iii) give rise to a population of interstellar and 
intergalactic MBHs. 

Subject headings: black hole physics - cosmology; theory - galaxies: nuclei - stellar dynamics 



1. INTRODUCTION 

The first massive black holes likely formed at high red- 
shifts (z ^ 10) at the centers of low- mass (~ IO^Mq) 
dark matter concentrations. These black holes evolve into 
the population of bright quasars observed at z ^ 6 and 
eventually grow into the supermassive remnants that are 
ubiquitous at the centers of galaxies in the nearby uni- 
verse. In popular cold dark matter (CDM) cosmogonies, 
dark matter halos and their associated galaxies undergo 
many mergers as mass is assembled from high redshift 
to the present. The merging - driven by dynamical fric- 
tion against the dark matter - of two comparable-mass 
halo+MBH systems will drag in the satellite hole towards 
the center of the more massive progenitor, leading to the 
formation of a bound MBH binary with separation of ^ 
pc. If stellar dynamical and/or gas processes drive the 
binary sufficiently close 0.01 pc), gravitational radia- 
tion will eventually dominate angular momentum and en- 
ergy losses and cause the two MBHs to coalesce. Such 
catastrophic events are one of the primary motivations for 
low-frequency, space-based, gravitational wave detectors 
such as the planned Laser Interferometer Space Antenna 
(LISA). For unequal mass pairs, gravitational waves also 
remove net linear momentum from the binary and im- 
part a "kick" velocity to the center of mass of the sys- 
tem. The dominant recoil effect arises from the interfer- 
ence between the mass-quadrupole and mass-octupole or 
current-quadrupole contributions (Peres 1962; Bekenstein 
1973). To date, the outcome of this "gravitational rocket" 
remains uncertain, as fully general relativistic numerical 
computations of radiation reaction effects during the coa- 
lescence of two Kerr holes are not available. For sufficiently 
asymmetric configurations, recoil velocities may exceed a 
few hundred kilometers per second and lead to a signifi- 
cant displacement of the MBH from the center of its host 
galaxy. In the shallow potential wells of small-mass ha- 



los at high redshifts, recoil velocities may be so large in 
the late stage of black hole-black hole coalescence to reach 
galactic escape velocities (Madau et al. 2004). If it is not 
ejected from the host altogether (giving origin to a popu- 
lation of intergalactic MBHs), the hole will return to the 
nucleus via dynamical friction. 

In this Letter we address the dynamics of recoiling holes 
in galaxy cores, discuss the implications of coalescence- 
induced kicks for the demography of MBHs, and consider 
the prospects for directly detecting the observational sig- 
natures of gravitational-wave recoil. While preparing this 
work for submission, we learned of an independent study 
by Merritt et al. (2004) of the consequences of the gravi- 
tational rocket. 

2. DYNAMICS OF RECOILING BLACK HOLES 

Gravitation radiation recoil is a strong field effect that 
depends on the lack of symmetry in the binary system. 
The lighter hole in a quasi-circular inspiraling orbit moves 
faster than the heavier one, and its gravitational radiation 
is more "forward beamed" . This gives a net momentum 
ejection in the direction of motion of the lighter mass, 
and the binary recoils in the opposite direction (Wiseman 
1992). According to quasi-Newtonian calculations (Fitch- 
ett 1983), at the transition from the inspiral to "plunge" 
phase, the center of mass of a compact binary of total mass 
M = mi -1-7712 (with the convention mi < TO2) recoils with 
a velocity 



vcM = 1480 kms" 



/max \ ^isco 



(1) 



= q2(l-q)/(l+g)5 

0.0179 for q = 1/2.6, 



where the function f{q = TO1/TO2) 
reaches a maximum value /max = 
and risco is the radius of the innermost stable circular or- 
bit, which moves inwards for a comparable-mass system 
relative to its test-mass limit (e.g. Buonanno & Damour 
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2000). By symmetry, the recoil vanishes for equal-mass 
nonrotating holes. Fitchett & Detweiler (1984) extended 
Fitchett's work to perturbation theory and estimated kick 
velocities > 100 kms~^. More recent perturbation theory 
calculations by Favata, Hughes, & Holz (2004) find that 
the recoil velocity can readily reach « 100 — 200 kms"-*^, 
but is unlikely to exceed « 500 kms~^. Numerical rel- 
ativistic calculations of radiation recoil from highly dis- 
torted Schwarzschild holes yields maximum kick velocities 
in excess of 400 kms"-'^ (Brandt & Anninos 1999). Because 
of the strong dependence of the rocket effect on the radius 
of the final prior-to-plunging orbit, larger kicks are ex- 
pected for prograde inspiral into rapidly rotating MBHs. 
Moreover, in the case of Kerr holes, recoil is significant 
even for holes of equal mass (Favata ct al. 2004) and can 
be directed out of the plane of the orbit (Redmount & 
Rees 1989). 

What is the dynamics of a recoiling MBH in the grav- 
itational potential of its host galaxy? As the MBH (with 
mass A^bh) recoils from the nucleus, stars boimd to the 
hole are displaced with it, so that the total ejected mass 
is Mtot = Mbh + Afcusp) where the bound stellar cusp has 



a mass Mr, 



Mbh for VCM ~ f, while Mcusp rapidly 



decreases for vqm ^ cr. As a first approximation let us 
assume that the hole -|- stellar cusp is on a radial orbit 
in a spherical potential (e.g., an early- type galaxy or the 
bulge of a spiral). Its orbit is then governed by 



GM{r) 



A-kG^ InApMt, 



ert(x) -j= 



(2) 



where x = v/^/2a, M(r) describes the mass profile of the 
galaxy, p(r) is the density profile of stars with ID veloc- 
ity dispersion cr, and the second term represents dynami- 
cal friction against the stellar background (e.g. Binney & 
Tremaine 1987; hereafter BT). We approximate the stellar 
density as an isothermal sphere, in which case the dynam- 
ical friction time scales with radius as tuF oc , while the 
'residence' time for a radial orbit scales as t k r/v (x r. 
Thus most of the decay of the MBH's orbit by friction 
takes place at small radii in the galactic nucleus, where 
it is reasonable to assume that the gravitational poten- 
tial is dominated by stars. Stars within the gravitational 
sphere of influence of the hole, i?BH ~ GMbh/c^, are 
bound to it and do not contribute to dynamical friction. 
We therefore truncate the stellar density inside a core ra- 
dius « i?BH, so that p{r) = (T^/[27rG(r^ 4- R'^yi)] ^'^^ 
M(r) = 2cr2r[l-(i?BHA)arctg(r/i?BH)]/G. Flattening of 
the; inner stellar density profile is also produced physically 
during the decay of the initial MBH binary, when dynam- 
ical friction and three-body interactions transfer energy 
from the binary to stars in the nucleus of the galaxy (e.g., 
Milosavljevic & Merritt 2001). 

The Coulomb logarithm In A in equation (2) is approxi- 
mately given by ln[6max/&min], where 6max and 6min are the 
maximum and minimum impact parameters for stars that 
contribute to dynamical friction (BT); femin is typically 
taken to be ^ i?BH (BT; Maoz 1993). In the case of a uni- 
form medium, 6,„ax is comparable to the size of the system 
(galaxy), so that 6max ^ ^min- For our problem, however, 
the decay of the hole's orbit occurs in the inner nuclear re- 



gions where the stellar density peaks, i.e. at radius ^ -Rbh, 
so that 6,„ax ^ -Rbh ~ '^min and InA ^ 1. Maoz (1993) 
has given a more careful derivation of dynamical friction 
in an inhomogeneous medium; the Coulomb logarithm is 
replaced (see his eq. 4.4) by In A jj^^^'' p{r)dr /(for), 
where i?max is the size of the stellar system and po is the 
central stellar density. For our problem this integral is 
^ 1, confirming the above argument that In A w 1 is ap- 
propriate. 



(T— 75 km s ^ 
M„„ = 3xl0'= M, 



v^y=400 km/s 




Fig. 1. - The response of a 3 x IO^Mq MBH to a given recoil ve- 
locity VQM in an isothermal potential with dispersion a = 75 km s~^; 
radial orbits are assumed. As explained in the text, the curves de- 
pend primarily on the dimensionless number vcm/c and so can be 
scaled to other parameters. 

Figure 1 shows the orbit of a 3 x 10^ Mq 
MBH moving with initial recoil velocity of vcm = 
200, 300, and 400 kms~^ in an isothermal sphere of disper- 
sion a = 75 km (reasonable for a Milky Way-type host) 
and radius i?max ^ ^bh- The escape speed from such a 
sphere is z;esc(0) = v'2|<^|(r = 0) = 2a ln^^\Rme../RBH)- 
We fix the hole mass using the Mbh ^ relation of 
Tremaine et al. (2002): Mbh = 1.5 x IO^ctIoo ^'-^0^ so that 
^BH ~ 15(t|oo where o"200 is the velocity dispersion in 
units of 200 kms~^. As Figure 1 shows, the timescale for 
the recoiling hole to return to the center of its host galaxy 
depends very sensitively on the magnitude of the kick it 
receives, ranging from ^ 10^ years for vcm = 200 kms~^ 
to more than 10^ years for vqm = 400 kms^^. Because 
the friction takes place primarily at small radii, the decay 
time is sensitive to the inner stellar density. For example, 
if we increase the core raditis of the isothermal sphere 
by a factor of 2, the decay times increase by a factor of 
w 3. Given uncertainties in the recoil velocity, this uncer- 
tainty in the decay time will not significantly modify our 
conclusions. 

The numerical results in Figure 1 can be understood an- 
alytically as follows. In the absence of dynamical friction, 
energy is conserved, and an ejected MBH will reach its 
apocenter at 



Ra ~ -Rbh exp {vcM/2a 



(3) 
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where the exponential dependence arises because the po- 
tential is logarithmic. During each passage through the 
galactic nucleus, however, the hole loses a fraction / of 
its orbital energy because of dynamical friction; this frac- 
tion can be estimated by comparing the DF timescale at 
~ i?BH to the transit time across this region, ^ -Rbh/''^cm7 
which yields / ^ (o"/?;cm)*- It takes ^ /^^ orbits through 
the nucleus for the MBH to return to the center of the 
galaxy, and so the decay time can be roughly estimated 
by tdccay ~ f~^Ra./vcM- Sincc both / and i?a depend pri- 
marily on wcm/c, this estimate shows that the decay time 
of the hole back to the center is determined primarily by 
vcm/ct for radial orbits. Thus, the results in Figure 1 can 
be scaled to other parameters, e.g. the decay times in Fig- 
ure 1 apply reasonably accurately to cr = SOkms"""^, so 
long as vcM is rescaled to vcm = 266, 200, and ISSkms"-^ 
(from top to bottom).^ 

The above analysis assumes that the MBH remains on a 
purely radial orbit. It is, however, likely that the recoiling 
hole will gain angular momentum upon leaving the nuc;leus 
(e.g. if the galactic potential is triaxial): this will modify 
the timescale for the orbital decay. Detailed estimates of 
the decay time in this case are uncertain. We can roughly 
bracket the effect by considering the MBH on a circular 
orbit with an energy (1/2)Mbhi^cM' ^'^''^ isothermal po- 
tential, the radius of such an orbit is « Raj \fe, where Ra 
is the apocenter radius for a radial orbit (shown in Figure 
1 and estimated in eq. [3]). The decay time is then given 
by the standard dynamical friction timescale (BT) 



10^'Vrs ( Ra 
In A 



1.5 kpc / V75 kms ^ 



(4) 



where we have used the Mbh — cr relation to eliminate 
the hole mass in favor of a. Note that the decay time is 
now determined by dynamical friction at large radii, rather 
than small radii, and so In A ^ 10 is probably more appro- 
priate. For vcyi « a the decay time predicted by equation 
(4) is comparable to that for radial orbits (a few crossing 
times; see Fig. 1) while it can be significantly longer for 
large vcm/o' (large apocenter distances). 



3. 



MBH BINARIES IN HIERARCHICAL CLUSTERING 
COSMOLOGIES 



In hierarchical clustering scenarios, MBH-MBH binaries 
are likely products of galaxy major mergers only. When 
two halo-l-MBH systems of (total) mass M and M, < M 
merge, the 'satellite' (less massive) progenitor will sink to 
the center of the more massive pre-existing system by dy- 
namical friction against the dark matter; for an isothermal 
sphere the Chandrasekhar dynamical friction timescale is 



top = 1-65 



1 



1 + MJM 



e 



(5) 



(Laccy & Cole 1993), where Avir is the density contrast at 
virialization, H is the Hubble parameter, and the term Q 
contains the dependence of this timescale on the orbital pa- 
rameters. After including the increase in the orbital decay 
timescale due to tidal stripping of the satellite (Colpi et al. 
1999), it is possible to show that satellites will merge with 



time only in the case of major mergers, Mg/M > 0.3 — 0.5. 
In minor mergers tidal stripping may leave the satellite 
MBH wandering in the halo, too far from the center of the 
remnant for the formation of a black hole binary. 

Major mergers are frequent at early times, so a signifi- 
cant number of binary MBH systems is expected to form 
then. It is still unclear whether halo major mergers nec- 
essarily lead to the coalescence of their MBHs, or whether 
the binaries "stall" before the backreaction from gravita- 
tional waves becomes important (Begclman, Blandford, & 
Rees 1980). A number of plausible mechanisms that may 
help avoid such stalling have been suggested in the recent 
literature (e.g., Gould & Rix 2000; Zhao, Haehnelt, & Rees 
2001; Armitage & Natarajan 2002; Yu 2002; Chatterjee, 
Hernquist, & Loeb 2003). 
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Fig. 2. — Mean number of major mergers (with mass ratios > 0.5) 
experienced per unit redshift hy halos with masses > 10^" M©- 
Solid line: progenitors of a Mq = 10^^ Mq halo at z = 0. Dashed 
line: same for Mq = 3 x 10-*-^ Mq. Dot-dashed line: same for 
Mo = 10^"^ Mq. The rates have been computed assuming the cur- 



rently favoured ACDM world model with fl^ 
h = 0.7, f^i, = 0.045, 0-8 = 0.93, and n = 1. 



0.3, = 0.7, 



Figure 2 shows the mean number of major mergers per 
unit redshift bin experienced by all halos with masses 
> 10^° M0 that are progenitors of a 2; = parent halo of 
mass Mo = 10^^ 3 x 10^^^ ^nd lO^^ Mg. We have tracked 
backwards in time the merger history of parent halos with 
a Monte Carlo algorithm based on the extended Press- 
Schechter formalism (Volontcri, Haardt, & Madau 2003). 
For the most massive parent halos this quantity peaks in 
the redshift range 2-3, the epoch when the observed space 
density of optically-selected quasar also reaches a maxi- 
mum (KaufFmann & Haehnelt 2000). Hydrodynamic sim- 
ulations of major mergers have shown that a significant 
fraction of the gas in interacting galaxies falls to the cen- 
ter of the merged system (Mihos & Hernquist 1994): the 
cold gas may be eventually driven into the very inner re- 
gions, fueling an accretion episode and the growth of the 
nuclear MBH. We will discuss below the possibility that 
gravitational-wave recoil may thus give rise to off-nuclear 
active galactic nuclei (AGN). 



the central galaxy on timescales shorter than the Hubble 

^More accurately, at fixed vcm/o', there remains a weak dependence on a via tdecay oc -Ra/^'CM oc -Rbh/'!'cm a. 
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4. IMPLICATIONS 

For galaxies with cr ;^ 50 — 75 kms~^, typical kick ve- 
locities of a few hundred km s"'^ are sufficient to unbind 
the hole or displace it sufficiently from the nucleus that 
the decay time due to dynamical friction is comparable 
to the Hubble time (Fig. 1). This implies that MBHs 
with masses Mbh ^ IO^Mq (using the Mbh ^ c relation) 
may be comparably rare in their late-type spiral or dwarf- 
galaxy hosts. Interestingly, there are very few observa- 
tional constraints on the presence of MBHs in such galax- 
ies, though Filippenko & Ho (2003) argued for a ~ IO^Mq 
MBH for the Seyfert galaxy in the late type (bulgeless) 
spiral NGC 4395, and Earth et al. (2004) reached a sim- 
ilar conclusion for the dwarf Seyfert 1 Galaxy POX 52. 
It is important to stress that even if galaxies do not cur- 
rently harbor a central MBH, they may have done so in 
the past. Thus the absence of MBHs in shallow poten- 
tial wells woiild not necessarily imply inefficient MBH for- 
mation (e.g., Haiman, Madau, & Loeb 1999), but could 
instead be due to recoil during MBH coalescence. Note 
also that the low-mass black holes that are preferentially 
affected by gravitational recoil are also those that are ex- 
pected to dominate the LISA gravitational wave signal 
from MBH-MBH coalescence. If ejection of such MBHs is 
common, this may decrease the number of sources detected 
by LISA. 

Galaxy mergers are a leading mechanism for supplying 
fuel to MBHs (e.g., Barnes & Hernquist 1991), and so a 
natural implication of gravitational recoil is the possibil- 
ity of off-nuclear quasar activity.^ Since the lifetime of 
merger-driven activity is of order the Salpeter timescale 
ts = 4.5 X lO'' yr, off-nuclear AGN are most likely to be 
found in relatively small potential wells. This is because, 
for (say) a lO^M© MBH in a cr = 200 kms'^ galaxy, 
the decay timescale due to dynamical friction is 10^ yr 
<C ts even for a kick velocity of ucm = 300 kms^^, and the 
displacement from the nucleus is quite small, ^ Rbh « 10 
pc. By contrast a decay timescale comparable to ~ ts 
is plausible for MBHs with Mbh ~ 10^ - lO^M© and 
i>CM ~ 200 kms^^. Since the merger rate for progeni- 
tor halos hosting such MBHs peaks at z « 2 — 3 (Fig. 
2), we predict that a significant fraction of moderate- to 
high-redshift, low mass AGN could be off-nuclear. Most 
currently detected AGN at high rcdshift are the rare 
~ 10* — lO^^M© holes hosted by massive halos, for which 
recoil is probably not important. At early times, however, 
theoretical models of the evolution of the quasar popula- 
tion in hierarchical structure formation scenarios predict 
a large number of fainter, Iow-Mbh sources (e.g., Haiman 
& Loeb 1998; Haiman et al. 1999; Volonteri et al. 2003). 

The effects of gravitational recoil may be particularly 



prominent in radio observations of AGN. Merritt & Ekers 
(2002) argued that MBH coalescence might imprint itself 
on the morphology of radio galaxies by changing the spin 
axis of the MBH and thus the direction of the radio jet. 
Gravitational recoil may have a comparably important ef- 
fect by displacing the radio-loud AGN from the nucleus of 
the galax;y. This could manifest itself as a flat spectrum 
radio core displaced from the optical nucleus of the galaxy, 
along with a jet or lobe still symmetric about the nucleus. 
Radio observations also have the obvious advantages that 
contamination from the host galaxy is less important (rel- 
ative to optical quasars) and VLBI interferometry allows 
precise localization of the radio emission. 

Gravitational recoil may also be relevant for under- 
standing the origin of off-nuclear ultraluminous X-ray 
sources in nearby galaxies (see, e.g., Colbert & Mushotzky 
1999). One interpretation of such sources is that they are 
intermediate-mass black holes (IMBHs) with M ~ 10^ Mq 
accreting near the Eddington limit. Miller & Hamilton 
(2002) suggest that such intermediate-mass holes could 
form in globular clusters by repeated black hole mergers. 
The gravitational rocket will, however, probably prevent 
substantial growth via mergers in the shallow potential 
well of a globular cluster, since even mergers with mass 
ratios as small as « 0.1 can lead to kick velocities in ex- 
cess of the escape velocity of the cluster. It should be 
noted that an IMBH could still form via stellar (rather 
than compact object) mergers during the core-collapse of 
young star clusters (e.g., Giirkan et al. 2004). 

Finally, we draw attention to the fact that the gravita- 
tional rocket does not necessarily limit the ability of MBHs 
to grow via gas accretion from rare less massive seeds, such 
as IMBHs produced by the collapse of Population III stars 
at 2; ~ 20 (Madau & Rees 2001). This is because seed holes 
that are as rare as (say) the 3.5a peaks of the primordial 
density field will evolve largely in isolation, as the merging 
of two (mini)halos both hosting a black hole is a rare event 
at these very early epochs. A significant number of MBH 
binary systems may form only later, when the fraction of 
halos hosting MBHs is larger. By then the typical host 
will be further down the merger hierarchy (more massive) 
and the effect of radiation recoil less disruptive (Madau et 
al. 2004). 
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